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ABSTRACT 


This thesis investigates the effectiveness of a search plan developed 
Бу В. 0. Koopman in a submarine versus Submarine search situation. Two 
computer simulation models allow probability of target detection as a 
function of sonar range to be used as a measure of effectiveness. The 
Koopman search plan is analyzed and a family of alternate search plans 
are developed. The choice of a particular alternate search plan is 
dependent on the parameters of the problem. These parameters are target 
speed, searcher speed, time late to the search area and total time avail- 
able to conduct the search. By use of the computer programs a search plan 
can be chosen so as to maximize the probability of target detection ata 


particular sonar range for each combination of input parameters. 
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І. INTRODUCTION 


During World War II a great deal of analytical and statistical work 
was done in the field of military operations. In particular the Navy 
benefited greatly from the work of the Operations Research Group (ORG) 
headed by Bernard 0. Koopman. The Group's research has been compiled into 
a single volume entitled Search and Screening (OEG Report No. 56). This 
volume explores the application of the scientific method to tactical prob- 
lems of naval warfare. It is intended to serve as a theoretical framework 
and foundation for more immediate practical studies. Philip M. Morse in 
the forward to Search and Screening (OEG ка No. 56) stated: 

"Although the tactical doctrines presented apply to 

instruments, weapons and conditions prevailing during 

World War II, it is believed that the methods and system- 

matic processes of analysis which led to the doctrines 

have wide application--not only to submarine warfare but 

to many other military and civilian problems." 

The analysis performed and the tactical applications devised by the 
Operations Research Group have withstood the test of time very well. 
Although extensive progress has been made in the methods of operations 
research since World War II, many of the tactical applications derived 
by the Operations Research Group are still considered the best available 
by current practitioners of the art. 

In particular the methods developed by Dr. J.M. Dobbie for constructing 
expanding square searches for targets in transit have been widely applied. 


However, some of these applications may be inappropriate. One use they 





have been applied to is that of a submerged submarine searching for 
another submerged submarine. In this particular situation the use of 
the search plan developed by Koopman and Dobbie is questionable. This 
thesis will investigate the appropriatness of the Koopman search plan in 


the submarine versus submarine search situation. 





Il. DISCUSSION OF THE KOOPMAN SEARCH PLAN 


A. BASIC ASSUMPTIONS OF THE KOOPMAN SEARCH PLAN 

In planning a search the nature of the target is usually known. The 
general position of the target is usually a random variable. Unless a 
fairly accurate estimate of its motion can be made, the search plan will 
have to be designed so as to be effective against a target having any one 
of many different sorts of motion. This plan assumes that the target speed 
is known exactly and that the target's course is a constant. When the 
object of search has nad its approximate location disclosed at a certain 
time, the searcher has the problem of disposing its subsequent searching 
effort (which is always limited) in such a manner as to maximize its chance 
of detecting the target. This is subject of course to the searcher's 
limitations and the practical limitations of navigation. 

It is assumed that information regarding the point of fix is received 
by the searcher. The time of this position information is also assumed 
known. It is further assumed that the searcher is airborne and has a con- 
siderable speed advantage over the target. 

The information concerning the target's position is not exact. Only 
a probability distribution of target positions is actually given. This 
distribution will have its greatest density at the point of fix and fall 
continuously to zero at a distance. The distribution will be assumed to 
be circular normal in form. Therefore, the probability that the target at 
the time of fix will be in a certain small region dx dy at the point (x,y) 


a distance r from the origin (point of fix) is given by the following: 





2 
-ү5/2 
f(x,y) = f(r} = 一 е ۴/26 





216 2 dxdy 


where 6 is the standard deviation of the circular normal 
distribution. 

It is assumed that the target's speed is known. The target's 
course is unknown but is assumed to be distributed uniformly between 
zero and 360 degrees. The target location distribution at a parti- 
cular time, T hours after the initial fix at a distance r from the 
origin (point of fix), can be expressed by the following probability 


density: 


Stak T 
L, >< — ) 





1 -(+r? 5 54217) /262) 


f(xysT) = 2C 
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The derivation of these relationships are shown in chapter two of 
Koopman [Ref. 1] where St is the estimated target speed. In this 
equation 1- is the ordinary Besel function of zeroth order value for 
pure imaginary arguments of the function. 

Figure I is a graph of 2r&<f(y;T ) this function plotted for 
various values of t as a function of r. It is obvious from figure 
one that when StT is greater than three sigma (36 ) the distribution 
has its maximum at r = StT. 

Because of this a search plan can be constructed for large values 
of time which have elapsed from the time of fix to the time of initia- 
tion of the search. This time interval is called time late. 

Theoretically the searcher should head for a position at a distance 


r where r = StT from the point of fix. This will place him on the peak 





of the target location distribution which is circular in form. The search 
pattern should then be constructed so as to maintain the searcher's posi- 
tion at the point of maximum density of the target location distribution. 
This would result in a search plan that resembled an equiangular or 
logarithmic spiral. There are practical considerations which make this 
infeasible, namely navigation problems and the undesirability of constant 
helming. Therefore, a search plan was constructed which approximates the 


theoretically optimal plan. 


8. CONSTRUCTION OF THE KOOPMAN SEARCH PLAN 

The theoretically optimal search plan was approximated by an expanding 
square search plan. Each circuit or square consists of four legs with 90 
degree turns. The shape of the squares are shown in figure 2. The leg 


lengths are E Эр” Las etc., and the corresponding distances of the 


2 
search legs from the point of fix are Г» Тоз Газ etc. The time required 
for the searcher to go from point A to point B is (ry + r2)/S a where So 

is the searcher's speed. The time required for the target to move from a 


distance r] from the point of fix to a distance r, from the point of fix 


2 
1S (ro - r4)/St where St is the target speed. The searcher will just keep 


up with the target if these times are equal. By this reasoning 


Үр = mr, 
Үҙ = mr, 
rq = mr, 
where 
о St 
ста 
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[+ we were to make ГБ = ТОД» the fifth leq would duplicate the position 

of the first leg in space relative to the target. That is to say that the 
searcher would be on the highest density point of the target location 
distribution. However, having traveled the first circuit (four legs) so 
as to be at this point of maximum density and not having detected the 
target has imparted additional information to the searcher. The first 
circuit being unsuccessfui has altered the target location distribution. 
The target location distribution after one circuit is no longer unimodal. 
It is now a bimodal distribution with a depression where the peak of the 
initial target distribution function had been. This is caused by the fact 
that one circuit has been completed without detecting the target. There- 
fore, the second circuit should be flown either de or outside the 
original peak. This will place the second circuit on one of the modal 
points of the new distribution. If the second circuit is constructed out- 
side the peak, the third circuit should be constructed inside the peak, $0 
as to be on the other modal value of the distribution. 

In order to accomplish this objective the fifth leg was moved a dis- 
tance S outside the position of the first leg in space relative to the 
target. This distance S is called the sweep spacing. A complete analy- 
tical development of an appropriate value for S is given in Ref. 1, chapter 
two and chapter seven. The final results of this development obtain the 


following form for the sweep spacing such that 


0.75 Ее 


Е = М2 1п2 


5 


мһеге 


= 
1 


sweep width 


2 x Sonar Range 


© = the standard deviation of the circular 
normal distribution. 
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Therefore, r- is determined so that 





5 
05" Пак Брег 5 
So | St 

so that 

r_ = mr4 + А 
where 

А = So х 5 

50 - St 


A first approximation to г" 15 St x T. However, the approximation to the 
spiral by straight search legs requires that г" be slightly less that 

this value. Therefore, a reasonable first approximation to Г] 15 obtained 
by setting rs 0.9 x St x T. Since any change in target course will 
reduce the outward component of velocity and since the second circuit is 


to be outside, the first г] was reduced further and taken so that 


гу * 0.8 x St x T 


Continuing we obtain the following: 


Го = mr. 
r3 = mr, 
rq 7 mra 
ГБ = ПГ +А 
rg = тг, 
r, = mre 
ra = mr7 
rg = mrg - ZA 





"19 = ША? 


TO 
ٹ٦‎ 
r13 = mrs» + 3 A etc. 


Using the relationship between the lengths of the legs and their dis- 
tances from the point of fix the following equations were obtained for 
the various search leg values: 


| „алш; 


1 1 

L, = ша + г. 
L3 = mL, 

ч = Ш. +А 
Г, = Щи 

Le = mle +A 
L; = ML. 

Lo = mL, - 2A 
Lo = ML, 
Lo = mio - ZA 
use 


| = е + 3A etc. 


This search plan was developed for large values of T where it was 
true that the maximum point of the distribution was located at a dis- 
tance r = St T. The question arises as to what is the smallest value 


of T for which this search plan may be used for all values of T larger 
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than this lower bound without any essential decrease in the probability 
of detection for a given amount of search effort. By observing figure 1 
it is seen that when T is less than &/St the distribution resembles the 
original circular normal distribution. For T larger than 26 / St the 
distribution has its maximum at approximately r = St x T and moves out- 
ward at St the speed of the target. However, there is a very rapid change 
in the distribution as T increases from €/St to 26 /St. This change of 
the distribution from a stationary distribution to a distribution moving 
at the speed of the target makes the decision on a lower bound for T most 
difficult. Koopman and his associates, by an empirical method of testing 
the search plan against various values of T, concluded that a reasonable 
lower bound on T was @/St. Therefore, for values of T26/St this search 
plan was considered optimal by Koopman and Dobbie. For values of T<6/St 


the plan for large values of T is sufficient and may be used. 


C. ANALYSIS OF THE KOOPMAN SEARCH PLAN 

It is necessary to evaluate the assumptions used in constructing the 
Koopman search plan when applied to its use by a submerged submarine 
searching for another submerged submarine. The assumptions of constant 
detectability over time and detectability being independent of target 
aspect are acceptable for the submarine application. However, the searcher 
obviously is "| атгрогпе пог does he possess a considerable speed advan- 
tage over the target. The pertinent question is this: Does the inapplica- 
ШІ. this assumption invalidate the search plan as devised by Koopman? 

In order to analyze this problem, a computer simulation model was 
designed. Various combinations of target speed, time late and searcher 


Speed were used to construct the Koopman search plan. By observing figure 





two, it is obvious that to be effective against a target whose course is 
assumed to be uniformly distributed between zero and 360 degrees, a 
minimum of one circuit or four search legs should be completed. The 
computer simulation constructed four legs of the search plan for each 
combination of target speed, searcher speed and time late. The program 
gave the time required to conduct one circuit of the search plan. In 
addition, it gives the probability of target detection as a function of 
actual sonar range as opposed to the estimated sonar range and resulting 
Sweep width used in constructing the search plan. In constructing the 
Simulation model, an option was included to allow the test of another of 
the basic assumptions of the Koopman search plan, namely, that of constant 
target speed. The program has two versions of target speed incorporated 
into it. In the first version the actual target speed is set equal to the 
estimated target speed used in constructing the search plan. The program 
can be changed so that the actual target speed is a random variable. In 
this case, the target speed is drawn from a uniform distribution. The mean 
of the distribution is the estimated target speed, and the end points of 
the interval are estimated target speed plus and minus two knots. This 
allows the evaluation of the search plan, when the assumption of constant 
target speed is not met, and an analysis of the resulting probabilities of 
target detection. Table I includes a tabulation of the amount of search 
time necessary to complete four legs of the Koopman search plan for various 
combinations of target speed, search speed and times late. 

As previously stated, the amount of search time available is always 
limited. From Table I it is apparent that the amount of time necessary 
to complete one circuit, herein considered the minimum against a target 
whose course is uniformly distributed between zero and 360 degrees, 15 


excessively large. The search time is particularly long for the most 


15 





probable cases that would be encountered in submarine versus submarine 
operations. 

Tables I and II allow a comparison of the probability of target 
detection for the case where target speed is deterministic as opposed to 
the case where it is a random variable. 

Because of the excessive amount of time required to complete the mini- 
mum number of legs of the Koopman search plan, it was decided to construct 
an alternate search plan that would not be susceptible to excessively long 


search times. 


16 











III. ALTERNATE SEARCH PLAN 


A. BASIC ASSUMPTIONS OF THE ALTERNATE SEARCH PLAN 

The Koopman search plan is obviously not universally applicable to 
Submarine versus submarine search due to the excessive search time 
required to complete even a minimum number of search legs. The cause of 
this failure is obvious. By examining the multiplicative factor m where 
5, + St 
ете 
the excessive time problem can be identified. 

For the range of target speeds and searcher speed combinations used 
in compiling Tables I and II, the factor m ranges from a high value of 13 
to a low value of 1.8. In the original design of the search plan, the 
searcher was airborne. Searcher speed was generally in the range of 140 
to 160 knots. The target, either a surfaced submarine or a surface ship, 
generally had a top speed of approximately 25 to 30 knots. These two 
facts taken together show that the multiplier m was severely limited in 
size. In the submarine versus submarine search problem, this is not the 
case. Ву examining the formulae for the various search legs it can be 


seen that the formulae reduce to a common form such that 


N 
L. = т гу + Cn 


where 2 (27: 


Therefore, for larger values of the factor m succeeding search legs 
become exponentially larger. This is the cause of the excessive times 


necessary to complete a minimum search plan. When target and searcher 
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speed are relatively close together and moderately large, this effect is 
even more pronounced. Therefore, an alternative search plan should not be 
limited in this manner. 

Another assumption implicit in the Koopman search plan is that numerous 
circuits will be completed. This fact was used explicitly in arriving at 
the form of the г. distance. The sweep spacing S was also determined based 
on the assumption that more than one circuit would be completed. Since the 
searcher in the submarine application does not enjoy the prerequisite speed 
advantage to complete many circuits, this assumption should not be used in 
constructing an alternative search plan. 

Intuitively then it seems that in an alternative search plan the multi- 
plicative factor would be relatively "small" so as to allow the completion 
of the minimum search plan in a reasonable amount of time. Also the r; 
distance would be determined based on the assumption that few circuits 
could be completed. This is due to the limited amount of Search effort 
available. In many cases tactical considerations will probably limit this 


time to that necessary to complete one circuit. 


B. CONSTRUCTION OF AN ALTERNATE SEARCH PLAN 

In order to test an alternative search plan against the Koopman search 
plan, a computer simulation model was developed. To allow a fair comparison 
of search plans, the search effort, that is search time, had to be equal for 
both plans. This time was set arbitrarily at four days (96 hours) for com- 
parison purposes. The target speed was eight knots, the search speed 12 
Knots, and the time late was ІІІ at four hours. The number of search 
legs of the Koopman search plan that could be completed in this time was 


then determined. The number of search legs for the alternate search plan 





was likewise determined. The probability of detection as a function of 
actual sonar range was then calculated for both search plans and compared. 
In determining the best alternate search plan, it was hypothesized that 
the construction would be as shown in figure 3. 
The distance г, would have the general form such that 
r] = (factor) x (time late) x (estimated target speed) 


where the factor value was as yet undetermined. The search leg lengths 


would be obtained as follows: 


L = (multiplier) x r 
n 1 


where the multiplier value was as yet undetermined. In determining compari- 
son Search plans, various combinations of factor values and multiplier values 
were used. The factor value ranged from 0.8 to 1.5 in 0.1 increments. The 
multiplier values used were 2.0, 3.0, 4.5 and 6.75. In all these values 
generated thirty-two different alternate search plans. 

The computer program output allowed a direct comparison of each alter- 
native search plan with the Koopman plan. The measure of comparison was 
the probability of detection as a function of actual sonar range. Although 
comparison values were available for all sonar ranges from zero to 60 nauti- 
cal miles, the comparison was made at a sonar range of ten nautical miles. 
This is the sonar range estimate used in constructing the Koopman search 
plan values of Tables I and II. 

The probability of target detection for different combinations of 
factor and multiplier values are shown in figure four. In this case, the 
actual target speed was equal to the estimated target speed. The Koopman 
search plan yielded a probability of target detection of .0870 under these 


values of target speed, searcher speed, time late and time available to 
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search. An alternate search plan using a factor value of 0.9 and a 
multiplier value of 2.0 produced the highest probability of target detec- 
tion (.1490). This is considerably better than the results using the 
Koopman search plan. 
However, this is not a fair comparison for the Koopman search plan. 
Koopman used the approximation for the г, offset distance as follows: 
г. = COX SL XL 


His previous approximation had been 


г. = 0.9 x St x TL 


The distance had been reduced because he felt that the target would not 
actually be on a constant course. The comparison is not fair because in 
the simulation model, the targets course is a constant. A better com- 
parison would result if the Koopman search were computed using an rı offset 
of 
ry = 0.9 x St x TL 

based on the fact that the target course is a constant once its initial 
value is determined. When this is done, the probability of target detection 
at a sonar range of ten nautical miles is .1140 for the Koopman search plan. 
Therefore, using a better comparison, the best alternate search plan is 
still better than the Koopman search plan. 

The question arises as to whether or not this particular alternate 
search plan is best or even better than the Koopman search plan when 
target speed is a random variable as opposed to being a constant equal to 
the estimated target speed. The simulation was run again with the actual 


target speed a random variable. The Koopman search plan yielded the 
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following values for probability of target detection at a sonar range of 


ten nautical miles: 


И 


r 0.8 x St x TL .0750 


1 
rı =0.9x5StxTL .0990 
The combination of factor and multiplier values that maximized the prob- 
ability of target detection at a range of ten nautical miles was the same, 
that 1s a factor value of 0.9 and a multiplier value of 2.0. They produced 
a probability of target detection equal to .1060. 

In theory what should now be done is to run the simulation for all 
interesting values of search speed, target speed and time late. This would 
be done with a set limit on the amount of search effort available. For 
eacn combination of target and search speed, the simulation would have to 
be run twice--once for the case where the actual target speed equals the 
estimated target speed and another time for the case where the target speed 
was a random variable. Each run would produce a different combination of 
factor value and multiplier value that maximized the probability of target 
detection. These values might be the same for the case where the target 
speed 15 deterministic and the case where it is a random variable or they 
might be different. One combination of factor and multiplier values might 
be ШЕПТІ for different target speed, search speed and time late 
combinations. This does not necessarily have to be the case. 

However, if this is not done, what is the applicability of the 
alternate search plan that was found best for this particular combination 
of searcher and target speeds and time lates with other parameter combina- 
tions? This hypothesis was tested by using the search leg values that 


this plan produces in the first computer simulation mentioned. The amount 
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of time necessary to complete four legs of the search plan for various 
combinations of search speed, target speed, and time late are tabulated 
in Table 1. 

Tables I and II list the probability of target detection for various 
combinations of search speed, target speed and time late values, using 
this particular alternate search plan. In the majority of the combinations 
the alternate search plan produces probabilities of target detection that 
are greater than the Koopman search plan. This is especially significant 
in light of the fact that the alternate search plan requires less time to 
complete one circuit than the Koopman search plan in all but three cases. 

Ideally, given a combination of search speed, target Speed, and time 
late, the alternate search plan that maximizes the probability of target 
detection for this combination should be determined. 

It is interesting to note a phenomenon exhibited by both the Koopman 
search plan and the alternate search plan. For some combinations of target 
speed, searcher speed and times late, the probability of target detection 
is higher when the target speed is a random variable than when the target 
speed 15 deterministic. This result occurs most often when target Speed 
and search speed are close together and the time late is large. 

It. is obvious that for some combinations of target initial position and 
a fixed speed, the target will be undetected by the searcher. However, if 
the target speed becomes a random variable, the searcher may be able to 
detect the target. Therefore, in these cases where the search plan has a 
low probability of target detection, randomness in target speed may in- 


crease the probability of target detection. 
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IV. RECOMMENDATIONS 


Since the amount of time required to conduct the minimum number of 
search legs using the Koopman search plan is excessive, it is recommended 
that an alternate search plan be used. The particular alternate to be 
used is dependent on the target speed, searcher speed and time late. The 
time available to conduct the search is also a deciding factor. The parti- 
cular search plan to use can be determined by using the m program 
described in appendix C. If this is not possible, the alternate search 
plan using a factor value of 0.9 and a multiplier value of 2.0 could be 
used. The resulting times to complete one circuit of the search plan and 
the probabilities of target detection are shown in Tables I and II for 


various combinations of searcher speed, target speed and time lates. 
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APPENDIX A 
TABLE I 


Constants for all Comparisons 
One Circuit (Four Legs) 
Actual Sonar Range - 10 nmi 
Estimated Sonar Range - 10 nmi 
Standard Deviation of Circular 
Normal Distribution - 2 nmi 
Actual Target Speed - Estimated 

Target Speed 


Koopman Search Alternate Search 
Plan Plan 


Time Searcher Target Search Prob. of Search Prob. of 
Late Speed Speed Time Detection Time Detection 
(HRS) (КТ5) (КТ5) (HRS) (HRS) 





3 8 4 160.61 | .4026 40.50 . 4586 
3 11 3 40.19 .4153 29.45 .4940 
СЕ 14 4 906 23.14 . 4150 
6 8 + 320.21 .1576 81.00 .1960 
6 11 1 79.81 .1880 22-11 .2010 
6 14 4 37.0 39272063 46.29 11383 
12 8 4 639.41 .0230 162.00 . 0783 
12 11 4 159.04 .0213 117.82 .0633 
12 14 4 ОЗ 02176 92257 .0600 
3 11 8 3416.89 0 68.91 .1670 
3 14 8 364.47 .1460 46.29 .1983 
6 11 8 6832.44 .0216 117.82 .0576 
6 1 8 HOC .0286 92757 .1006 
12 11 8 13653-53 .0000 235.64 . 0096 
12 14 8 1,455.87 0 185.14 20225 
3 14 12 64,026.36 .0773 69.43 МІРОЮ 
6 14 12 128,050.69 .0000 138.86 .0656 
12 14 22275650996 5.0000 277.71 10123 
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TABLE II 


Constants for all Comparisons 
One Circuit (Four Legs) 
Actual Sonar Range - 10 пт 
Estimated Sonar Range - 10 nmi 
Standard Deviation of Circular 
Normal Distribution - 2 nmi 
Actual Target Speed - Uniform 

Random Variable 


Alternate Search 


Koopman Search 
Plan 


6 


12 
12 


12 


Plan 
Searcher Target Search Prob. of Searc Prope OT 
Speed Speed Time Detection 1те Detection 
(KTS (KTS) (HRS) (HRS 
8 4 160.61 .3616 40.50 .4033 
11 4 40.19 .4000 29.45 .4040 
14 4 19.20 .4313 23.14 .4033 
8 4 320.21 .1113 81.00 . 1266 
11 . 4 79.81 .1190 58.91 . 1403 
14 4 37.99 . 1486 46.29 .1310 
8 4 639.41 „0340 162.00 10320 
11 4 159.04 .0376 117.82 .0356 
14 4 75.59 .0453 92.57 . 0373 
11 8 3416.89 „1290 58.91 . 1903 
14 8 364.47 „1373 46.29 „1786 
11 8: 6832.44 .0330 117.82 0723 
14 8 728.27 .0456 92.57 .0643 
1: 8 13,663.57 .0086 235.64 .0213 
14 8 1,455.87 .0166 185.14 .0193 
14 12 64,026.36 .0810 69.43 .1463 
14 12  128,050.69 .0190 138.86 .1516 
14 12 256,099.63 .0026 277.71 . 0186 
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APPENDIX В 


This appendix discusses the first of the two computer simulation 
models used. This simulation determines the amount of time necessary to 
complete one circuit (four legs) of the Koopman search plan. The output 
of the simulation is a graph of probability of detection as a function of 
actual sonar range and a listing of the sixty points plotted by the graph. 
The program was written in Fortran IV, level G and was run on an IBM 360 
digital Computer. ж 

The basic program organization is discussed next. The basic input 


parameters that can be changed are as follows: 


ALPHANUMERIC 

NAME PARAMETER VALUES UNITS 

NTRUN Number of targets run >] 
against a particular 
search plane 

RS Sonar Range 60>RS>0 NMI 

SIGMA Standard deviation of >0 NMI 
circular normal distribution 

50 Initial value of searcher 51 «50414 KTS 
speed 

STE Initial value of target 0«5Т«50 KTS 
speed estimate 

TL Initial value of time 0<Т1.<12.0 HRS 

` Тате 

ULCT Upper limit on target O<ULCT <21 rad 
course 

LLCT Lower limit on target 0411.СТ«21 rad 
course 
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For this program the input parameters mentioned above were initialized 


as follows: 


NTRUN = 3000 
RS = 10.0 
SIGMA = 2.0 
SO = 8.0 
SIE = МЕО 
TL = 4.0 
ULCT = 21 
ЦЕСТ = 0.0 


The program then uses these initial values to compute the parameters 
necessary for constructing the Koopman expending square search plan. 
The four search leg lengths are then computed. These distances and the 
г] distance are then used to establish the (x, y) co-ordinates for each 
Search leg end point. Using the search speed input, the time at which 
each of these co-ordinates is reached is determined as is the total search 
time necessary for the entire search plan. The headings entitled initial 
values for search legs one, two, three and four will be discussed after 
the determination of target motion is discussed. The target course is 
then obtained from a uniform distribution. The end points of this interval 
are те as previously mentioned. In this case the target speed was 
selected uniformly between values of zero and 360 degrees. The target 
speed is then obtained from a uniform distribution. The actual target 
speed can be set at the estimated target speed by using the following input 
cards: 

ULST = STE 

ELST = SIE 
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Or the actual target speed can be a random variable drawn from a uniform 
distribution whose end points are the estimated target speed plus or minus 
two knots. This is done by substituting the following cards in the input 
deck: 

ULST 


STE + 2.0 
LLST 


5ТЕ - 2.0 

The direction and distance of the uncertainty in target position from the 
origin (point of fix) is then generated from a circular normal distribution. 
The program then uses the randomly selected target course to determine the 
target's velocity components in the x and y directions. Using the random 
target position uncertainty selected from the circular normal distribution, 
the initial (x, y) coordinates of the target are then determined. 

The conceptual framework associated with the next section of the pro- 
gram is discussed in detail in this section. For each randomly selected 
target track there will be some time at which the target and searcher are 
at a minimum range. This range is called the closest point of approach 
or CPA range. 

It is useful to think of each search leg as being extended in space 
past the point where the turn to the next search leg is made. At any 
time on any search leg the searcher's position is a function of time and 
his ОР) That is to say: 

(ај + ant,a, + agt)‏ = لی لاہ و×) 
where‏ 


search leg initial x- coordinate 


> 
= 
II 


as = search leg Initial y- coordinate 


ap = searcher x- velocity component 


ад = searcher y- velocity component 
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In a similar manner the target's position can be expressed as follows; 


(х+ Y+) = (by T bot 4 ba T bt) 


where 
b. = target's initial x- coordinate 
b. = target's initial y- coordinate 
b = target's x- velocity component 


by = target's y- velocity component 
It also can be shown that the searcher/target range is obtained as 
follows: 
2 2 tie 
Range = (x -x,)° + (ys-yt) 
Now to minimize this range the derivative of the equation with respect 
to time is taken and then is set equal to zero. This procedure yields 
the time at which the CPA range occurs. This time is given by the 


following equation: 


esd (41-51) (42-02) + (аз-В3) (ад-бд) 
(42-85)2 + (ag-bg) 

If this time actually occurs while the searcher is traversing this search 
leg, the range at this time is a candidate for a minimum range. Ir the 
minimum range on this search leg does not occur on the interior of the 
search leg interval, then it must be an end point solution. The target/ 
searcher range, when the searcher is at the end points of each search leg, 
is then easily determined. 

The target's speed is a constant once selected. Therefore, the x,y 
velocity components are constants. The target's initial position is also 
constant once determined. On a particular search leg the searcher x,y 


velocity components are constants. One is identically zero, depending 
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on the orientation of the search leg. Therefore, for each search leg 
there are at most three candidates for the CPA range and at the least 
two. If the solution to the time equation falls within the time frame 
that the searcher is on this particular search leg, then there are three 
candidates for a minimum on this search leg. These three are the search 
leg end point ranges and the range at the time equation solution. If 

the time value that satisfies the equation occurs when the searcher is 
not on this particular search leg then there are only two candidates for 
a minimum range, the search leg end point ranges. Using this procedure 
all possible candidates for minima can be obtained for the entire search 
plan. Once the CPA range for the iteration has been determined, it is 
placed in a tally box. There are sixty tally boxes. Each tally box 
corresponds to an interval of one nautical mile. Once the program has 
determined the CPA range for each target, this range is placed in the 
tally box for the interval into which it falls. A11 CPA ranges in excess 
of sixty nautical miles were placed in the last tally box. The program 
tests at this point to determine if the required number of targets have 
been run against this particular search plan. If not, the program gener- 
ates another set of target parameters and repeats the above mentioned 
step. If the required number of runs has been obtained, the program sums 
all the tally box values and then determines the probability that the CPA 
range will fall in the interval of each of the individual tally boxes. 
This information allows the calculation of the probability that the CPA 
range will be less than or equal to a particular range. In essence this 
is the probability of target detection as a function of the actual sonar 
range. The output of the program is twofold. First, a print out of the 


probability of target detection as a function of sonar range is provided. 
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These values are listed for every nautical mile between one and sixty 
nautical miles. Secondly, the probability of target detection as a 
function of sonar range is plotted graphically for values of sonar range 
from zero to sixty nautical miles. The program then varies the searcher 
speed, target speed and time late to obtain all the combinations shown in 
Tables I and II. 

The probability of target detection as a function of actual sonar 
range is in fact an empirical cumulative distribution function. The num- 
ber of iterations performed of different target tracks (3000) ensures that 
the empirical cumulative distribution function will be within + .025 units 
of the population cumulative distribution function. This is with a prob- 
ability of .95, that is to say % = .05. The exact number of runs needed 
was arrived at using Kolmogorov's statistic. These values are tabulated 
in Table A21-b of Ref. 2. 

The pseudorandom number generator used was tested for randomness using 
a separate computer program not described in the thesis. The statistical 
tests for randomness used were the frequency test, serial test, lagged 
product test, runs up and down test, and runs above and below the mean 
test. The pseudorandom number test passes all of these tests adequately, 
thereby justifying its use. 

It should be noted that by replacing the г. distance calculation and 
the search leg length computations for the Koopman plan by those necessary 
for the alternate search plan, the same output can be obtained for the 


alternate search plan. The following card must be removed: 


SM = (SO + STE)/(SO - STE) 
r= GUC SIE TL 
SLI = SM * г. 


зі 





512 


ӘМЕ ۶۲ r. 
513 = SM * SL2 
SL4 = SM * SL3 +A 
Once these cards have been removed the following cards have to be inserted 


in their places: 


SM = 2.0 

гу 40.9 * STE * TL 
SL] = SM * г, 

SL2 = SM ** 2 * r. 
1 MS Te 
SL4 = SM ** 4 * r. 


The program will then determine the values depicted by Tabies I and II 


for the alternate search plan that was determined to be best. 
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APPENDIX C 


This appendix describes the second computer simulation used in the 


analysis and comparison of search plans. The primary purpose of this 


simulation is to compare alternate search plans with the Koopman search 


plan. Therefore, the goal of this simulation is to compare search plans 


which require the same amount of search effort, that is search time. 


With this goal in mind, the input variables were chosen as follows: 


ALPHANUMERIC 
NAME PARAMETER VALUES 
STE Estimated target speed 0<5ТЕ<50 
50 Searcher speed STE<SO 
TL. Time late >0 
TLimit Amount of search time »0 
available 
NTRUNS Number of target runs »0 
against a particular 
search plan 
ULCT Upper limit target course O<ULCT<2N 
LCT Lower limit target course O<LLCT<2H 
SIGMA Standard deviation of >0 
circular normal distribution 
RS Sonar Range >0 


The program was run with the following initial values for 


STE = 8.0 
50 = 12.0 
TL = 4 
NTRUNS = 3000 
ULCT 2т 
LLCT = 0.0 


83 


UNITS 
KTS 
KTS 
HRS 
HRS 


RADIANS 
RADIANS 
NMI 


NMI 


these parameters: 





2.0 


516МА 
RS = 10.0 
The program input considered next is the actual target speed input. 
There are two options available. First the following two cards may be 
used: 


ULST 


STE 

LLST = STE 
When used, these cards result in the actual target speed being a constant 
equal to the estimated target speed. The second option involves inserting 


the following cards into the deck: 


ULST = STE + 2.0 

LLST = STE - 2.0 
This results in the actual target speed for each iteration being drawn 
from a uniform distribution. The mean of the distribution is the estimated 
target speed while the end points are the estimated target speed plus and 
minus two Knots. The program then computes the parameters necessary for 
the calculation of the Koopman search plan. The (x,y) coordinate values 
for each of these search leg end points is then calculated. This having 
been done, the values for twenty search legs were computed. The values 
necessary to determine the searcher's initial position and velocity com- 
ہش‎ or each of the search legs is computed. The program then determines 
the time that the searcher will reach each of the search leg end points. 
The procedure stops the first instance in which this time exceeds the limit 
of the search time available. By so doing, the number of search legs or 
parts of search legs that can be completed in the search time available 
has been determined. The target's course is then selected from a uniform 


distribution. The limits of this distribution are variable and have been 
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previously mentioned. The target's position uncertainty is then selected 
from a circular normal distribution. The next target parameter determined 
15 target speed. This speed comes from a uriiform distribution as mentioned 
before. 

The target's initial (x,y) coordinates are then determined. The tar- 
get's velocity components a the x,y direction are calculated next. The 
candidates for the minimum target/searcher range are then calculated using 
the procedure fully explained in appendix A. From these a CPA range for 
this target iteration is determined. This range is placed in a tally box 
and the procedure is repeated until the required number of target iterations 
has been performed. When this has occurred, the probability of target 
detection as a function of sonar range is calculated. 

The input parameters for the alternate search plan are now listed. 

They are as follows: 
OMULT 
FACTOR 


2.0 


0.8 


The factor value is used in determining the г" distance as follows: 


ШЕ- FACTOR * SHES TE 
and the quantity OMULT is used in determining the search leg lengths 
where 


ж n 
L. = OMULT (гу) 


The value of twenty-four search leg's lengths for each FACTOR and OMULT 
combination are then determined. The (x,y) coordinates values for each 
search leg are calculated. The parameters necessary for the target/ 
Searcher minimum range calculation are determined at this point. The 
number of full and partial search legs that the target can complete in 


the allotted search time is then determined. The target course, speed 
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and position uncertainty values are then determined. The closest point 

of approach is determined for this iteration. This CPA range is placed 

іп а tally box. The program tests to see if the desired number of targets 

have been run for this search plan. If they have not, it generates tar- 

gets until the required number has been reached. For each target it com- 

putes the CPA range and places it in a tally box. The program then computes 

the values of probability of target detection as a function of sonar range. 
These values are then printed for the Koopman search plan and the parti- 

cular alternate search plan under consideration. The program then plots 

the graph of probability of detection as a function of sonar range. The 

values for the Koopman search plan and the particular alternate search plan 

under consideration are then plotted on the same graph. This is done to 

facilitate comparison. The factor value used in the г, determination and 

then the multiplier value used in the search leg length calculation are 

then varied. The entire alternate search is then recalculated. The 

required number of targets are run against this search plan and the prob- 

ability of detection as a function of Sonar range is calculated. Once 

these values are calculated they are printed out and plotted for easy 


comparison with the Koopman search plan. 
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COMPUTER OUTPUT 


A sample of the computer output is included in this section for each 
of the programs used. The output for the first program consists of a 
listing of the probabilities of target detection for actual sonar ranges 
for a limited number of target speeds, searcher Speeds and times late 
combinations. After each listing, the same information is presented in a 
graphical format. The x-axis represents actual sonar range in nautical 
miles, and the y-axis is the اس‎ of target detection. In using 
the listing, a computer design feature made the following convention 
necessary: 


Probability of Target Detection = CDF(x+1) 
for Sonar Range of x nmi 


The output of the second computer program consists of comparative 
listings and comparative graphs for the input parameters specified in the 
program explanation of Appendix C. The listings and graphical results 
are given for the Koopman search plan and the alternate search plans. 

The alternate search plans used a multiplier value of 2.0 and factor values 
of 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, and 1.5. The same convention is used 
for the listings. The axes of the graphs are also the same. In the graphi- 
cal representation of the results, the Koopman search plan values are 
denoted by “." and the alternate search plan values are denoted by "+". 

In the case where the values of the Koopman search plan and the alternate 
Search plan under consideration are the same, only the value for the 
alternate search plan (+) is printed. The exact probability of target 


detection at a particular sonar range can be obtained from the listings. 
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МОГ ОУЧУСО со OM SMO MM OO OLN OO NNT me O UO mH OD TVDOOLTLUIINOOONGOONTNNTOOOOODOOODO 
VOON TM O MUI O Mi MON $ U UO N یہ ہہ‎ QUI DO ہہ‎ HUN O DUNA UA LA ATA سر‎ ОО ОоОосососоОосооОосСосооОоОоО بت‎ 

СОЕ ہے‎ NN NANO а зет FT FIONN O OT TT ODO SOOO O O i O со 
со о © ө ө ф e с с с с с с с с с с о с с с с с с 6 с с с с с с с с с с с со с с с с с с с с со с с е ® ес с с с с е өе ө 
ОоООООдсООООООСООоосооооосоооооооосоооосооооооооооооооон ہے ہے ہں ہے مس ہے ہے ہب ہے ہے‎ 
l H H ٢ H H H H H HI lH H H H H H I H W H H H H H H l H H H H H ٢ ٢ ا ا‎ U ا‎ A H I H H H H H ٦ H HH H H H H l H и 
می +ہ یہ بے‎ О یہ ‌بہ ن ح‎ ۰+ М со О INIA YT تا‎ MA со СО mm NIN) LAS = OQA O ян A Tuy O> O ÇA O —rN n T uy افو‎ OO یں‎ = 
یہ یہ یہ یہ یہ یہ یہ یہ یہ ہے ٭ى- ہہ ہے ہے ہے ہے نے ہے ہے‎ ہل)٦۱‎ MN OOO TIN N TIN YF f Tq ^ dq ما ۱۳ما ۱۱ای خ بج ی سپ‎ ~ 


UL UL Lis Li» Lis Li LL, Li» Lie Lis UL UL UL. LL. UL LL. UL. UL Lis Li» Lis Lis Li Li» Li. Li» Li» Lis Li. Li» Li. Ú U U U. UL. Ú. U UL. U UL UL. U. U U. UL UL U. UL UL UL U. U. Li. Uis Li. Li» Li» Lis Li. Lis 
БЈЕЈЕТЕЈАТЕТЕТА ГЕТЕ ТАТЕТЕТАТЕТЕТЕТЕТЕТА ГЕТО ТЕТЕТЕТЕАТЕТЕТЕТЕГЕТЕТЕЈА ТЕЛЕ ГЕТЕ ГЕТА ГЕТЕ ГЕТЕ ТА ГЕТЕ Те ГЕТЕ ТЕ ТЕ ТЕЛЕ ГЕТЕ 15] 
СООООООООООООООО ООО ОООООО ООО ОДО ОДО ОДО ОДО ДОО О ن ں ن ن یں‎ о Ооо 
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OM MOO OR CORE EER CEO Or Om O mm O mm ars Сге ГУ РУ ГУ О (ГУ О посопоси- сссоссосососзо 
OI O N O QO MM O O 0 СО О ОЗО یی ص۹ بت‎ )۱ са соса 40 са О 42 ога ооо по ото ОД Ооо о Ооо Ооо Ооо О 

A Y (N SD DDO DU FO SO SMA OND ODO 09 OO (MM y MN Ч O O DN чо счча AC CAO O Оо О О О О О О 

O O" (MW O FO QUID IO OO mM CUP™ UN N CO UE CO N ON SJ O OY (NIYA DO A Y ےئ ںب‎ = CO CO UA ОХОО ООУ СО Со یاں یا‎ CO بں‎ AG 
Шатты атанады аге гг" со مں دی دی‎ ODONTOL OOOO A کاپ‎ 
• е е е е ә е ә е в о ә 9 ә е о е 9 е ә е е ә е о ә о е ә о е о о о ә о ? о ә е е е е е е 0 6 ә 6 0 0 0 0 0 в а ө в е е 
ТОСЕ О И ООО ОТ ж ОЮ Боо O O کا ہے ہے ہے ہے ہہ ہے نہ ہے ہے رک ن تان کات نا نات کا یق ری کا‎ 
Hl H H jJ H H H H H H H H RE lo H H H H H H H H H ll H ов ов H H H l H H HI H H H H Io I H H HH I H H H H H I H H H H H H H H 
A NAO +f iX QO 0 PO ںیئ +۷ لہ‎ S OO TF O —cN cy -T'UYAQ - O QY O — یہ‎ ۱٣7ی‎ Pa 09 DBO SAIC LAY DAA رں‎ UNO — rN rO ELA ГО СУ О sd 
اسم‎ і үний pu п о тай ај аф С) СМ) СОЈ Су Су СЧ Су СУ Су ССО СО СО СУ СО сос т хот SES OS LALA LALA LA О О 
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шы ندانادناناتیا)نانانانادناناناتاناتانائتاناناتانادناتاتاناناناتاتاناناتن نادناتاناتاداتنادناتنانادادنانادننانادنانانادانانادا) انا ت٦ا دا‎ 
ҮӚХхххх ن ی ی ی یں ی یں یں ی ی‎ Y OA IO OU OO ن ی ن ی ن ن ی ی ی ی ی یں ی یں یں ن یں ن ی ن ن ی ی ن ی یں یں ی یں یں نے‎ DIO VOW 
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ССО ГУ ОСО ОФ О Се М түе) С СГУ М СО О СОГУ ГУ СО СО Се 62 Г Со СО СО СО МСО М СО О со соса ГУ К ОЗ MMO NMV OV MNO 
HHNIIDMOODOMIIANAMMOOGOOIMIDIMOMOIMMOMIOMLUMOMMMIIODODMMOINDDODOMUIO 

GAT Or- Qo‏ ہمان OO f UN O mm P~ OO CN ONE EAL Qe LALALA UO S O Y MO O О СУГЛОО Г TÁ NO TW TONTO‏ بے 

COCO ں ۳ہ یہ ن‎ DOWN UI A ی۹۳ یہ ہے یی ںو نہ‎ Y 09 DO AULA ге со ہے بی یی‎ NN O O) УУГУ ГУ СО сол ОО а О 
OIOOVOIOYDO OO OO O I ۷ہ) لہ یہ لہ یہ یہ یہ یہ ہے ہے ہے ہے ہے بے ہے ہپ‎ AI AAN کب‎ I ISIS А АГАР АГАР НУО 
60.0 0 e е е о о о 9 в می‎ с с о = 0 0 06 0 0 0 0 6 0 6 ә е ә ә е е е е а ее ә Ө е е ә е е ә ә ә е ә ә е ә е е е ә P о в в е 
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H ll H HI HH H H HN H H HH H H H H HH H H H HM HH HH H H lH HII H HH HN H H HMHH H HI HH HMHI JM H HH HH HI 
СЧС УО М СОСО Чс PLA O = OO OA یہی- ن‎ ( О М со С о СЧС О М ООСО ЧОЧ СОМ ооо со ихо г< оо соч 
لح لہ بہ ہے ہے ہے ہے ہے ہے ہے بے ہے ہے‎ N NINA III IN NINO ТЭ T F T YT ٣۳ ۱ماما گا ما گا گا١کا گا تاب‎ O O 


ш.Ш.Ш.Ш.Ш.Ш.Ч.Ш.Ч.ШМ.Ч.Ч.Ч.М.М.Ш.Ч.ЧЦ.Ц.Ц.Ц.Ч.Ш.Ч.Ш.Ц.Ц.М.ЦШ.М.Ц.Ш.Ц.ШМ.ЦМ.Ц.Ц.ЦМ.Ц.Ц.Ш.ШМ.Ш.Ц.Ц.Ш.Ш.Ш.ШМ.Ц.Ц.Ц.Ш.Ц.М.Ц.Ц.М.ЧМ.Ц.Ч. 
(JC ود‎ ыс ыс کت جب 656.263 65 ٤ب ت اب 653 626 62 6565 6.3 6265 62 6.265 65 653 62 62 626.262 62 62 26.263 62626265 6262 62 62 62 62 6262 6.26262 جح بی‎ 
ل ن ن ن ن ن ن ن ن ن ن ن ل نل نل ل ن ن ن ن ن ن ن ن ن ل ن ن ن ن ن ن ن ن ن ن ن ن ن ن ن ن ن یں ل ن ن یں ن ی ںی ن ن ن یں یں یں ن ن یں‎ 
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87.58 


OM Orn OU Ore № сп OO O OOM O ОМ М m m O n Or r n rr MIMMOrTMMMMDONTOOmr О 
MOIDOMMOODDDOM DVO VUMOLOY DOU MOMOUMOMMMOVUMMODMMMMOOMVOOVO 

O)‏ بی ہے ری mMN (NO O M N O O m UN O mU‏ ن (O Y UOT 3 СЪХГ СУ СО СО PPh OMOEA AT ON eI‏ ںہ سہ FF‏ یىی COO‏ ١کالہ‏ ہم .حم 

Фс чао оччого О чеч риа ар N mH ION ОР ОС NNT FF О ГУ UO N O HN У ЧУЧО ГУ ООЛ ОСЧ г O 

2999999 09990090444 ۷اا یہ یہ لہ یہ یہ‎ NTI EQS TE TTT TT TAINAN INN OD YOO СЭ 
е 9 еее е ее еее еее ее ее ее еее е еее еее о о о 6 е е е “ е е е е е е е ев е е “ е جج‎ е е вз еее в 
а ое осо ооо особо бо соо о бобообо сова 
UU U UN O U DUN U HH H HH DU YH U UN U ۷۱۰ UH U U U ٢ ا‎ HN NH H N H H HW HNH WH H H IH H NN H H HNHH HHN HH H HNH H 
ہی ۹۴ یہ بے‎ со ن ہ وہ ح مہہ یہ بے ں ی‎ АС MN ГУ ао О О mH NN ЦУ У ن ی وو‎ С AND О TO HNN TINO ON Om 
یہ لہ یہ نہ لہ یہ یہ یہ یہ یہ ہے ہے ہے ہے ہے نے ا ہے ہے ہے‎ ٦١۱١٠۲۱۱۳۹ ۱٢۹.۰ N rr f oq dq dq گا اہ ہی۴‎ ۱۱ NO 


ц.ЧЦ.Ц.М.М.Ц.Ц.М.м.Ц.Ч.Ч.Ш.Ц.Ц.Ш.М.Ч.м.М.Ц.М.М.Ц.Ц.М.Ц.Ц.ЧЦ.Ч.М.М.Ч.Ц.Ц.М.М.Ц.Ч.Ц.Ц.М.Ц.М.Ц.Ц.М.Ц.Ч.Ц.ЧМ.М.Ц.М.М.М.ЦШ.Ц.М.Ч.й Lhu 
ذا بت ںے بت تباب بب اف‎ NOOO OOOO OO OOOO QUO QV UNO NOV OO HOODOO 36336244 303303 andino 
эссэ نا ںان‎ ә نا ںان نال لاہںنا ںا ںا ں نب ں ل ل ں نں‎ ООСС ОДО СО О ООСС ОЈ 
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OOOOOOOOOOOO OO VOOOVOVOOOOVVOVOOVOVODVOVOOOVOVOOOVOIYVOVOOVOVOYOD 

Y MP O OIE NIN ON MSO LIB 00 O N) O CY mU CY OO OD یی‎ A N O J" mm ری‎ ОО СМ нні NY A АР ХР У Oun O) LAMA I OO 
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نف ی یں ں ب ہے ں 9وی 1 ںی یں ں ہیں ای ان ںب‎ о ں ں و ویو ن0 0 نں ںی ب ن‎ C OCO ہے وہ00 0ب و 0ں نں ب نں ن0 ٗن ۱ ں0۱ 0 ں ہو مر‎ o o ہن‎ 
l H H H H H H H H UO O ٢٢ ا٢۱‎ H H HH H W H Ó H H H H H H H UO H NH H H H H H H ۱٢ H W ۱٢۲ H H HH M H H NO H H H H Ó W HIH H и 
MAN MO SPUN OP OD ORONO) SUD OM OD ON MISE LN OM OD OMAN MO МОГ O TO mM NN $ ЧА ЗО Г оо ХО ч NN ЧУО ПГ UNO rd 
(۱۳۰؟۲() لہا لہ لہ لہ لہ یہ لہ) پہ یہ لی ہے ہے ہے ہس ہے ہے نے بن بے نے‎ )۴( ۲٦ ۱٦ )۲۱۳۴٦( ۴٦ р .گا گا ا ما ۱۳۱ ۱گا ۴ا۱گا ۱گا ہی ۳ ٣پ عي -ب ٣پ يٻ‎ О 
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O OO O O O Q OUO O O O یں‎ O OY OO OO O O OU O O O O YO NO O O O O O ی‎ NOOO ODO O O O VO DO DO OD VO WY 

пі ПО А ЧУ Ге РФ UA Л СОСЫ UA ДУМА С ГИ АЈ یہ یں ہے یہ‎ CA Y O(N ON mM OP NY" OU OD ں‎ DINO TM O UD NV O 
OOONT wd MON MUY DO + y یئ‎ DD یہ ہے‎ р ШҮ مم یہ بت یىی‎ ST 0 PM 00 MO O NRO SD Y Pe DA Ц УГ ON OO Pe ن ی ں‎ 

O O O OUO Q O O O O O =m mind mM NN NN NNN 1°10 9 C9 إ7(‎ © Y3 J” “پ‎ d٣ .ب‎ ہ٣‎ ХР АР ۱ا۴١ کا گا‎ ٦اک‎ LY ГМ ГАЧ УО Чу) М) GOTT OAHU 
ее ее е ее е ө е е 6 9 е е 6 ө е в ә е е о о е 6 ә ә е е ев е е е ә е е е е 6 е е 6 е е е е е е е е е е 6 е ә е ө е е в 
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یہہ‎ )۷٦( ٣ یہ بت یں رہ ۔٢ رک ا‎ ЈСО YT UO رہ‎ N Om ری نہ ریہ یں ٠ہ یں‎ N O Се СОЧ) ГУ ОО Ол OD ама ЧУО => OOO — (N (Q Np UA Г а) UNO) md 
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ہداب ы ыгы сзсз ы сзсз сзсзсзсзсо ш‏ ب بات بب ہ۱ تب 62 بن ت ыссы‏ جب ہج جت ہت 2262© 622 С3‏ 26262 دج تب جح ہ بب جب جح ت جب بنا تی 
ی ل ں ن ںین پب ایی اب ب ن یں ن ن ن ں بہا‌ٴہپراں یاہاں پ نپ ن ن ن ن ین ن رای نان با دژاںن0|اہان ذانب اب ن ںیاب ں تاہیاں باب ن ل ں 
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ОООООООООООООООДОООООООООоОооооооооооооооооооооооооооооооооооо=~ 
H HM H H ٢ ٢٢ H ll M H H H H H N HNH HÓ W ۱٢۰ H HI H H H HW H I H IH H H IH H H H H N ٣٢٢ ٢ ٢٢ ٢٢ ٢ ٣٢ ٢٢ ٢ ٢٢ ٣٢ ٢ ٢٢ا‎ HN H ٢ H 
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W H H Il Hl H H H H H H WRH H H HN H H NH N H H H H l MH H H H H ٢٢ H ٣۲۰٢٢ ٢٢ HH H N H H H H H H U MO H ON HOHH H H H H M H 
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РЦ У = s en e s G G G 0 £ o s G é 6.4 98 G QO 9 9 8 Ө 9 8 8 86 6 8 8 6 6 e 0 G ө а өө © э өе ө © е э ө GR GS G q ос 
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U- Ue U Ue Lt Li U LL LL Lh Lh, Li LL, LL LL LL LL, LL, LL, LL LL LL U LL LL LL LL LL L LL LL UL tL LL LLL U لا ما ا لا لا تا‎ ушу ا لا لا نا نا‎ 
ОСсспспсоспсопопосососоососзоспооссосососзосоосопбоооосос) ہب تب ہے نہ دا یجباجح دنا نے ہ یناب چنا‎ 
ОКО ОКО OW O ن‎ OO VO VO OO OO OO OO VO OO OO یں ن ن یں یں ن نن ن ںں ن ن ںں یں ں ن ں‎ 
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